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Nitrogen-vacancy centers in diamond, being a promising candidate for quantum information processing, may 
also be an ideal platform for simulating many-body physics. However, it is difficult to realize interactions 
between nitrogen-vacancy centers strong enough to form a macroscopically ordered phase under realistic tem¬ 
peratures. Here we propose a scheme to realize long-range ferromagnetic Ising interactions between distant 
nitrogen-vacancy centers by using a mechanical resonator as a medium. Since the critical temperature in the 
long-range Ising model is proportional to the number of spins, a ferromagnetic order can be formed at a temper¬ 
ature of tens of millikelvin for a sample with ~ 10“* nitrogen-vacancy centers. This method may provide a new 
platform for studying many-body physics using qubit systems. 

PACS numbers: 76.30Mi, 76.50.+g, 05.50.+q 


Introduction. The negatively charged nitrogen-vacancy 
(NV) centers in high-purity diamond have been considered 
as a promising candidate for solid state quantum informa¬ 
tion processing due to their long coherence time CHSl and 
high feasibility in initialization, control, and readout of their 
spin states 0. Simulation of many-body physics using NV 
centers, in analogue to cold atom physics 10151, has been 
proposed CHU- To realize phase transitions in NV cen¬ 
ter qubit systems under realistic temperatures, however, suf¬ 
ficiently strong interactions between NV centers located suf¬ 
ficiently close (<30 nm) is till highly challengingllTl [8l [TTII . A 
new opportunity is to use resonators as mediators llT HT^fTBl . 
which has potential of coupling NV centers at distance. 

In this letter, we propose to realize long-range coupling be¬ 
tween many separated NV centers via a mechanical resonator. 
A remarkable feature of the long-range interacting system is: 
the critical temperature for ferromagnetic phase transitions is 
proportional to the number of spins, so a ferromagnetic or¬ 
der could be formed at a temperature of tens of millikelvin 
(mK) for a sample with ~ 10^ NV centers, even though the 
mediated coupling between two NV centers is less than 200 
kHz(~ 10 yuK). 

Microscopic model. The essential idea can be understood 
by considering a prototype system shown in Fig.[2a). A me¬ 
chanical resonator with frequency cOr is attached with a mag¬ 
net. The NV sample is placed right under the magnet with a 
distance d. The oscillation of the mechanical resonator gener¬ 
ates a time dependent magnetic field that causes Zeeman shift 
of the spins in the sample. In contrast to the schemes of direct 
dipolar coupling between NV centers, our scheme does not re¬ 
quire short distance between NV centers. We use a diamond 
sample where the orientation axes of the NV color centers can 
lie along one of the four possible crystallographic axes in dia¬ 
monds, [111] (Fig.[2b)) which can be grown by chemical va¬ 
por deposition method nulla. The NV centers have a spin 
triplet ground state (5 = 1) with a large zero field splitting 
A = 2.87 GHz. We apply a static magnetic field Byvv along 


the [111] direction of the sample so that two Zeeman states 
of the electron spin qubits with magnetic quantum number 11) 
and |0) are nearly degenerate, 6 = (A-yeByw) ~ 0 (Fig.[TJc)). 

We make use of the two Zeeman states of the electron spin 
11) and |0) and define a pseudo-spin by cr^ = 11)(1| -10)(0| and 
the corresponding spin flip operators cr^ = |1)(0| and cr" = 
|0)(1|. In each NV center there is a '^^N nuclear spin which 
interacts with the on-site NV center spin. The Hamiltonian 
of the NV center near the degenerate point (d ^ 0) is //nv = 
2;[Aiv(/7)^ - jNl'jBNv + Alj ■ a-j -t- doj], where A^ = 5.1 
MHz is the '"*N nuclear spin quadrupole splitting 1151, JN is 
the gyromagnetic ratio of the nitrogen nuclear spin and A « 2 
MHz is the hyperfine coupling between the electronic spin and 

the *^N nuclear spin II51- 

The motion of the mechanical resonator is described by the 
Hamiltonian Hr - ojrb'^ b, with cjr as the fundamental vibra¬ 
tion mode of the resonator and b as the corresponding annihi¬ 
lation operator. For example a silicon nitride string resonator 
has dimensions (325 x 0.35 x 0.1 )pm with a>r = 27r x 1.0 MHz 
and Q- 1.3 x 10^ at room temperature |[T71 . 

The magnetic field felt by the electronic spin of the NV cen¬ 
ters can be approximated by a magnetic dipole IfTSlI . B{dj - 
x(t)) - B(dj) - Gmx{t) -I- o(x^), where dj is the distance be¬ 
tween the equilibrium position of the resonator and j-th NV 
center, Gm is the magnetic field gradient at the position of the 
NV centers and x is the amplitude of mechanical resonator os¬ 
cillation (x ~ 10“'^ m at temperature of mK). This magnetic 
field will induce a Zeeman shift to the NV spins with Hamil¬ 
tonian //j = ^nd rjj - jelGmlao, where jg is 

the gyromagnetic ratio for electron and ao - ylftllmtOr is the 
amplitude of zero point fluctuations for a resonator of mass m. 

Due to the distance distribution of the NV centers to the 
magnet, the coupling of the NV centers to the mechanical res¬ 
onator have a distribution. A magnetic tip with size of 100 
nm produces a magnetic gradient G,„ ~ 7.8 x 10® T/m at a 
distance 25 nm away from the tip |[T9l . In such case rflln 
could reach 200 kHz. For a sample with typical distance of 
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FIG. 2. (color online). Free energy landscapes of the long-range 
Ising model, (a). The free energy as a function of the spin polariza¬ 
tion along the .r direction at high temperature = NJlhgT = 0.5. 
(b). The same as (a) but at a lower temperature ^ = NJIkgT = 2.0. 


FIG. 1. (color online), (a). Schematic set up of a hybrid NV and 
mechanical resonator system. A magnet attached to the end of the 
mechanical resonator is positioned at a distance 25 nm above the NV 
centers, coupled to the electronic spin of the defect centers, (b). The 
NV defects in the samples. The defects have one preferred orien¬ 
tations in the sample, [111]. A static magnetic field is applied 
along the [111] direction, (c). Schematic energy levels of NV centers 
as a function of the magnetic field. The Sjvv is taken so that the states 
|1) and |0) of the NV centers are nearly degenerate. 


NV centers ~ 20 nm, the neighbor NV centers have a direct 
dipolar interaction with strength ~ 5 kHz (~ 0.25 yuK), much 
less than the critical temperature to be discussed later. For 
the sake of simplicity, in the following we assume the cou¬ 
pling between the mechanical resonator and NV centers to be 
uniform and neglect the direct dipolar coupling between NV 
centers. Without these assumptions, however, the results in 
this paper would only be quantitatively affected. The cou¬ 
pling constant rjlln ~ 200 kHz considerably exceeds both 
the electronic spin coherence time (milliseconds) of the NV 
centers and the intrinsic damping rate of the mechanical res¬ 
onator, y - Ur IQ of high-2 mechanical resonator Ifm . We 
neglect the interaction between the mechanical resonator and 
the ^"^N nuclear spins since they are three orders of magnitude 
smaller than that between the NV center and the mechanical 


resonator. 

Summarizing the three terms described above, the total mi¬ 
croscopic Hamiltonian including the spin qubits in the NV 
centers, nitrogen nuclear spins, the mechanical resonator and 
the interaction among them is, 


H = Urb^b -I- ^ [AyvC/p^ - jnBnvI'j + Alj ■ o-j\ 

J=i 

N 

+ ■ ( 1 ) 

M 

Effective long-range Ising model. To obtain the effective 
interactions between NV centers, we calculate the partition 
function of the hybrid system Z = Tr[e“^^]. We shall neglect 
the Zeeman energy of the '“^N nuclear spin since its strength 
(~ 0.3 MHz) is much smaller than the zero field splitting ~ 5 
MHz and the hyperfine interaction ~ 3 MHz. The hyperfine 
interaction between the electronic spin and the '^N nuclear 
spin term AI ■ cr can be separated into the diagonal term AFcr'‘ 
and the flip-flop term A{PcB’ + Ecr^). Let us first neglect the 
flip-flop term and the qubit splitting (d = 0). By partial trace 
of the phonon bath and the nuclear spins, the partition function 
can be factorised as 


Z = Tr[e-^"] ~ ■^phonon ^ ^ Tr[exp( 
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FIG. 3. (color online). The phase transition temperature of the NV- 
resonator hybrid system as a function of the scaled detuning d/Wr. 
We take the NV-mechanical resonator interaction strength 77 = 200 
kHz, eigen frequency of the mechanical resonator co,- = 1 MHz and 
the number of spin is = 10000 . 


where Zphonon is the partition function of an independent 
phonon bath, is the partition function of the nuclear spins 
and the effective Hamiltonian for the NV electron spins is 


‘ 4 —^ COr 
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with hifi) - ln[(l + 2 cosh(;6Aj))/3] being an effective 
magnetic field resulting from the on-site interaction between 
the NV center and the nuclear spin. The physics for the 
effective long-range interactions between the NV centers is 
that an NV center spin flips by virtually emitting a phonon, 
and then the phonon is virtually absorbed by another NV 
center spin far away with a spin flip. Therefore two NV 
spins at distance can flip-flop via virtual exchange of phonons, 
which leads to long-range cr^crj coupling. In the tempera¬ 
ture range we are interested in(~ 10 mK), the effective held, 
hij3) ~ 0.3//K, is negligible. Therefore the distant NV centers 
interact through an effective long-range ferromagnetic Ising 
interaction with the Hamiltonian 

^electron “ “ ^ 

i<j 


where the ferromagnetic coupling strength J - Irf'lur. If 
the flip-flop term of the hyperflne interaction is taken into ac¬ 
count, it would contribute an effective transverse held to the 
electronic spins, which, however, is ~ 0.3/7K and hence negli¬ 
gible. 

Ferromagnetic ordering. The long-range Ising model can 
be solved by mapping the ensemble onto to a single large spin. 
With the notation Sx - the Hamiltonian of the 

long-range Ising model can be written as 

//electron = ^ " f)' 

i<J 
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FIG. 4. (color online). Comparison of the finite size numerical re¬ 
sults and the mean field theory, (a). The order parameter, Mx, as a 
function of the reduced dimensionless temperature, kgT/NJ for dif¬ 
ferent sizes of the spin systems. The black-short-dashed line is for 
N = 100, the blue-long-dashed line is for N = 500 and the red-solid 
line is the mean field prediction, (b). The specific heat plot in the 
same style as in (a). 


The spin quantum number S takes values 0,1,2, • • • , N/2 for 
even N and 112,312,- ■■ , N/2 for odd N. Moreover, the spin 
degeneracy in each large spin subspace is /)[5] = - 

^ * where is the binomial coefficient. Then the elec¬ 
tron spin partition function 


N 

n=0 


In the large N limit ll20l 


Z 





( 6 ) 

(7) 


where ip = (1/2 - 1 - X)ln(l/2 + X) + (1/2 - Z)ln(l/2 - X)- 
2pX^ and p - pNJ. The integration in the partition function 
can be evaluated by the saddle point approximation method. 
For p < \, the saddle point appears at Z = 0, and for p > 
1, there are two symmetric saddle point located respectively 
within (-1/2,0) and (0,1/2). Therefore = 1 is the critical 
point for the ferromagnetic phase transition in the long-range 
Ising model. So the critical temperature for the long-range 
ferromagnetic Ising model, ksTc - NJ , is proportional to the 
number of spins and the coupling strength ||2T1 . 
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We present the scaled free energy of the long-range Ising 
model as a function of the spin magnetization = {S x) IN 
for different temperatures in Fig. One can see that at high 
temperature (Fig. [^a)) the free energy minimum is situated 
at Mx - 0 and a paramagnetic state is stable, while at low 
temperature (Fig. |^b)) the free energy minimum is doubly 
degenerated with Mx - 1 or Mx = -1, and a ferromagnetic 
state is preferred. 

For the hybrid system, the phase transition temperature can 
be controlled through the detuning from the degeneracy point 
between the NV spin states |1) and |0). If the NV qubits have 
a small splitting 5 = A - yBi^v, the resonator mediated inter¬ 
action between the NV centers becomes 


- 


2 2 
6 OJr + 


i)' 


The critical temperature becomes 

Nt]^ Nrf- 


ksTc - 


Ur - 6 


Ur 


( 8 ) 


(9) 


Here the magnetic field is tuned so that the NV spin states 
|1) and |0) are near degenerate. Hence the splitting 6 is much 
smaller than the critical temperature. In Fig. we plot the 
phase diagrams of the NV centers as a function of splitting 
for N - 10000. In a diamond sample with NV concentration 
~ 50 ppb, coupling between a resonator and such many NV 
centers is possible. A ferromagnetic order can be formed at 
the temperature of 50 mK. 

In Fig. Qa), we present the numerical and the mean field 
results of the order parameter as a function of reduced temper¬ 
ature, ksT/NJ. One can see that the order parameter vanishes 
at and above the critical temperature hgT/NJ = 1 and also the 
finite size numerical results approach to that of the mean field 
predictions reasonably well as the number of spins increases. 
The specihc heat shown in Fig. j^b) has a finite jump which 
signals a second order phase transitions occurs. 

Conclusion. In summary, we propose a scheme to real¬ 
ize ferromagnetic ordering of distant nitrogen-vacancy cen¬ 
ters by using a mechanical resonator to mediate long-range 
Ising-type interaction. The critical temperature for the fer¬ 
romagnetic phase transition in the long-range Ising model is 
proportional to the number of spins, so the ferromagnetic or¬ 
der could be formed at the temperature of tens of millikelvin 
for a sample with ten thousand nitrogen-vacancy centers. In 
addition, it may also be possible to use a superconducting res¬ 
onator Eaiia as a medium to realize the long-range ferro¬ 
magnetic coupling between the NV centers. Since the interac¬ 
tions between the NV centers mediated by a superconducting 
resonator is usually small compared to that by a mechanical 
resonator, high density NV samples are required to observe 
the magnetic ordering in the NV centers. 
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